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1. Introduction 

The intergalactic magnetic field (IGMF) is a weak, possibly primordial field that permeates 
the voids between galaxies [1]. One method to probe the strength of the IGMF is to use obser¬ 
vations of distant blazars with imaging atmospheric Cherenkov telescopes (lACTs). A fraction of 
the TeV-range y-rays emitted by blazars will undergo pair production upon interaction with 
the extragalactic background light (EBL). These secondary particles will then interact via inverse 
Compton scattering with low-energy cosmic microwave background (CMB) photons, up-scattering 
the photons to higher energies. To illustrate the energies involved, a ~10 TeV photon emitted by 
a blazar will be reprocessed to secondary cascade photons of ~100 GeV energy [2]. Thus, it is 
possible that if a blazar has a sufficiently hard spectrum and sufficiently high energy cut-off, cas¬ 
cade emission could be observed by an lACT arrays such as VERITAS. The currently operating 
lACTs have energy thresholds around 100 GeV, thus for an appreciable fraction of the emission 
above 100 GeV to be due to cascade emission, the source must have a cut-off at TeV energies. 

In the presence of a non-zero IGME, the cascade emission will be angularly broadened and 
time delayed as the deflection of the intermediate electrons and positrons by the magnetic field in¬ 
creases the path length between the source and the observer compared to unscattered emission. The 
magnitude of the angular broadening and time delay varies with the IGME strength and can be di¬ 
vided into three regimes. Eor 10^^^ G < B < 10^^ G, the electron-position pairs will be isotropized 
in the vicinity of the blazar, forming a pair halo detectable as an angular broadening compared to 
unscattered emission and the lACT’s point spread function (PSE) [3]. Eor 10^^^ G < B < 10^*^ G, 
the bulk of the pairs will not be isotropized. The cascade will be slightly time delayed as well 
as angularly broadened compared to unscattered emission [2]. Eor B < 10^*® G, the predicted an¬ 
gular broadening is too small to be disentangled from the PSE, at least for the angular resolution of 
the currently operating lACTs, and the cascade can only be detected via observation of a time de¬ 
lay compared to unscattered emission [4]. This study focuses on searching for angular broadening 
rather than time delays, and is thus insensitive to the third regime of IGME strengths. 

In addition to the dependence on the intrinsic source spectrum, the projected sensitivity to 
cascade emission depends on the source redshift, due to the evolution of the EBE with redshift. Eor 
very distant sources (z > 0.2), the cascade emission is too compact to be distinguishable from the 
beamed emission, whereas for nearby (z < 0.1) sources, the cascade emission is too broad to be 
distinguishable from the cosmic-ray background. 

Previous searches for extended emission around blazar have been performed by MAGIC (using 
Mrk 501 and Mrk 421), Fermi-LPW (using a large blazar sample), VERITAS (using Mrk 421) and 

H. E.S.S. (using lES 1101-232, lES 0229-t200 and PKS 2155-304) [5, 6, 7, 8]. 

I. 1 Source selection 

The sources used in this analysis were selected for optimal sensitivity to magnetically-broadened 
emission. As stated above, the predicted cascade fraction is greatest for hard-spectrum sources with 
high energy cut-offs. The spectral index F measured by Fermi-LKV in the GeV band is used to 
approximate the intrinsic index [9]. Sources with r<2 were selected. Note that in the absence 
of a Ferm/-EAT-measured spectral index for lES 0229-1-200, an intrinsic spectral index r=1.5 is 
assumed. The cut-offs in the intrinsic spectrum were taken from the joint GeV-TeV fits performed 
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z 

r 

Ecu,[TeV] 

^detect 

Mrk421 

0.031 

1.771±0.012 

0.692 

185.3 

Mrk501 

0.034 

1.738±0.027 

1.630 

94.8 

H 1426-t428 

0.129 

1.316±0.123 

0.570 

7.6 

VER J0521-t211 

0.108 

1.928±0.034 

10 

23.2 

lES 0229-t200 

0.14 

1.5 

10 

10.3 

lES 1218-t304 

0.182 

1.709±0.067 

1.010 

35.5 

PG 1553-tll3 

0.5 

1.665±0.022 

10 

46.0 


Table 1: Source list, with the redshift, the assumed intrinsic spectral index (based on the Fermi-LAIT mea¬ 
sured index), high energy cut-off, and detection significance in the dataset used for this analysis. 


in [10]. For sources not appearing in [10] or consistent with no spectral break, a cut-off energy of 
10 TeV assumed. This value was selected based on its position at the lower edge of the range of 
cut-off energies that could still yield a spectrum consistent with no spectral break. 

To minimize statistical uncertainties, strongly-detected sources were used, ideally those with 
> 10a significance. This requirement was softened for sources falling in the redshift range with the 
greatest expected sensitivity (z~ 0.1-0.2). Although sources at the optimal redshift were preferen¬ 
tially selected, sources at higher and lower redshifts were also included in the analysis. In the case 
of a detection of a broadening of the emission, this would allow tests of the predicted redshift de¬ 
pendence. The final source lisf, fogefher wifh fhe assumed infrinsic specfral index, energy cul-off, 
redshifl and defecfion significance in fhe dafa sample used is shown in Table 1. 

2. VERITAS observations 

VERITAS is an array of four 12 mefer lACTs [11] locafed af Fred Lawrence Whipple Obser¬ 
vatory in soufhern Arizona. VERITAS is sensifive fo very high energy /-rays in an energy range 
from 85 GeV fo greafer fhan 30 TeV. VERITAS has an energy resolufion of 15% and an angular 
resolution of 0.1° af 1 TeV. 

The data used in this analysis were collected between 2009 and 2012, prior to a major camera 
upgrade. As obtaining the best possible angular resolution is critical for this measurement, only 
runs with all four telescopes operating were used. Additionally, only runs with zenith angle of 
observations <30° were used, both because the angular resolution is better and because the energy 
threshold is lower than for large zenith angle observations. During strong flares, if is expecfed 
fhaf fhe unscaffered emission will be brighfer fhan fhe cascade emission, decreasing fhe overall 
sensifivify to cascade emission. Thus, periods of high source acfivify were removed from fhe 
dafasefs. 

3. Analysis 

All dafa were processed using fhe sfandard VERITAS calibrafion and shower parameferiza- 
fion mefhods. The /-hadron separafion was achieved using selecfion on mean-scaled Hillas pa- 
ramefers for all sources ofher fhan H 1426-1-428, lES 0229-1-200, and VER J0521-I-211. For fhese 
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Figure 1: Comparison between the angular profile of three sources and their simulated point sources. Good 
agreement is clear from the residual distributions. The results of a probability test are shown in Tab. 3 for 
all sources. 


three sources, which were not detected to a high significance with standard y-hadron separation, 
a boosted decision tree (BDT) analysis was used. The BDT analysis was not used for the other 
sources because the energy threshold is slightly higher for the BDT analysis, and preserving the 
lowest possible energy threshold was prioritized for sources already strongly detected with the 
standard y-hadron separation. 

To achieve a good angular resolution while retaining as many photon events as possible, events 
with images recorded in only two telescopes were discarded, while events with three and four 
images were retained. The angular resolution is additionally correlated with the distance between 
the shower core and the array center. The selection on this parameter was optimized to improve the 
angular resolution, resulting in an optimal value CoreDist <215m. 

To search for extended emission, the angular distribution of a source can be compared against 
a known point source whose angular extent is due only to the instrument PSF. The parameter 6^, 
defined as the angular distance between the air shower’s arrival direction and source’s estimated 
location, characterizes the angular profile. The VERITAS PSF varies with energy, the azimuthal 
and zenith angles of observations, and the night sky background level. Thus, it was necessary to 
simulate a point source for each blazar studied that matched the source characteristics and observ¬ 
ing conditions. The range and mean value of the night sky background level was matched between 
each source and its simulated point source. The simulations were weighted event-by-event to match 
their data counterparts’ energy and azimuthal angular distributions. It was not possible to exactly 
match the sources’ zenith angular distributions, as simulations were generated at discrete values 
for the zenith angle of observation. A zenith angle of 20° provided the best approximation of the 
zenith distribution in data, and was therefore used. A correction (described below) accounts for 
this simplification. 

The agreement between the 6^ distributions for the sources and their matching simulated point 
sources was assessed first based on the derived residual distributions and a probability test. 
Examples of several overlaid data/simulation distributions are shown in Figure 1. To produce a 
measure of the agreement that includes systematic uncertainties and the zenith angle correction. 
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Wdata i (Tstat 

^sim i ^stat i ^syst 

Mrk 421 

0.0483° ± 0.0004° 

0.0484° ± 0.0002° ± 0.0009° 

Mrk 501 

0.0503° ± 0.0008° 

0.0481° ± 0.0003° ± 0.0007° 

lES 0229-t200 

0.0429° ± 0.0075° 

0.0461° ± 0.0003° ± 0.0011° 

H 1426-t428 

0.0591° ± 0.0109° 

0.0547° ± 0.0003° ± 0.0014° 

VER J0521-t211 

0.0428° ± 0.0004° 

0.0451° ± 0.0002° ± 0.0012° 

lES 1218-t304 

0.0468° ± 0.0020° 

0.0507° ± 0.0003° ± 0.0007° 

PG 1553-tll3 

0.0558° ± 0.0021° 

0.0521° ±0.0002° ±0.0011° 


Table 2: Fitted width of the 9 ^ distribution for sources and their corresponding simulated point sources. 


the angular distributions were fit and the width parameters Wdata and Wsim were compared. 

3.1 Systematic uncertainties 

Although the dominant uncertainties in this analysis are statistical, the uncertainties intro¬ 
duced by mispointing and energy resolution on the simulated point sources were assessed. The 
uncertainty on the telescope pointing is ~25", which translates to an uncertainty of 0.0005° on the 
6^ distribution. 

The PSF depends on energy, thus the finite energy resolution of VERITAS introduces a further 
systematic uncertainty on the 6^ distribution. An average energy resolution of 20% was assumed. 
The propagated uncertainty varies depending on the steepness of the source spectrum, and was thus 
calculated for each source. The uncertainties vary between 0.0005° and 0.0014°. 

As mentioned previously, the simulated sources were generated assuming a zenith angle of 
observations of 20°, which only approximates the observed zenith distribution in data. This as¬ 
sumption was addressed by deriving a function PSF (Ze) for the PSF as a function of zenith angle 
using a large sample of Crab Nebula observations. A PSF width Wze{obs) was then derived based 
on the observed zenith distribution, and the difference {wze{obs) ~ wze=20“) was used to correct the 
simulated source’s fitted width. The size of the zenith corrections range from 0.0009° to 0.0023°. 

4. Results and Discussion 

The fitted widths derived from the source and simulated point source 6^ distributions are 
shown in Tab. 2. Although some marginal tension is present - Mrk 501 appears to be broader than 
a point source at the 1.95a level and lES 1218-1-304 appears to be naiTower at the 1.87a level - 
there is no significant discrepancy between the sources and their simulated point sources. 

Checks of the agreement using the 6^ histograms based on residual distributions and prob¬ 
ability tests also indicated good agreement between the sources and their simulated point sources. 
The p-values corresponding to the x^ probability tests are shown in Tab. 3. Only PG 1553-1-113 has 
a marginal p-value, which does not account for the zenith correction or systematic uncertainties. 

4.1 Model-dependent limits 

Model-dependent limits were set for assumed IGMF strengths of B = 10^'^ G, 10^^^ G and 
10^^^ G. As extended emission was not detected, magnetic fields in this range provide the most 
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p — value 


Mrk421 

0.19 

Mrk 501 

0.38 

lES 0229-t200 

0.30 

H 1426-t428 

0.95 

VER J0521-t211 

0.31 

lES 1218-t304 

0.52 

PG 1553-tll3 

0.003 


Table 3: Results of a comparison of the 9^ histograms for the sources and their simulated point sources with 
a probability test. 

straightforward interpretation. In the case of non-detection, it is not possible to rule out magnetic 
field strengths in the regime in which the pairs are isotropized and a pair halo forms without 
making strong assumptions about the source’s history of activity. 

As discussed in the introduction, the projected sensitivity to the cascade emission hinges heav¬ 
ily on the source’s intrinsic spectrum, and in particular to high energy cut-off. For VER J0521-I-211, 
lES 0229-1-200, and PG 1553-1-113, Ecut is taken to be 10 TeV and the fraction of the total emis¬ 
sion due to cascade emission is within VERITAS sensitivity (cascade fraction>5%). Eor the other 
sources, the predicted cascade fraction is too small to warrant setting model-dependent limits. 

Model predictions were obtained from a 3-dimensional semi-analytical cascade code [12], 
using the Gilmore EBE model [13] and assuming an IGME coherence length of 1 Mpc. The blazar 
jet’s Doppler factor was set to 10, and the viewing angle was set to 0°. The simulation took into 
account the source redshift and events were weighted to match the desired intrinsic spectrum. 

An ensemble of pseudo-experiments for the extended emission was generated by convolving 
the predicted extended emission 6^ distribution with the PSE distribution corresponding to the 
source of interest. The uncertainty on the point source width was incorporated by sampling a 
Gaussian distribution with a width reflecting the statistical and systematic errors on Wsim (hence 
the use of pseudo-experiments rather than a single convolution). A second ensemble of pseudo¬ 
experiments at each cascade fraction fc was generated. The appropriate number of y-ray events for 
each source was sampled according to [(1 —fc)x Direct emission H- /^x Extended emission]. 

Eig. 2 shows the fitted width of pseudo-experiments for lES 0229-1-200 versus fc, overlaid 
with the width measured in data. The uncertainties on the measured and simulated widths, added in 
quadrature, are considered in the extraction of upper limits on fc. Upper limits were derived for PG 
1553-1-113 and VER J0521-I-211 as well, however lES 0229-1-200 allows the strongest constraints 
on fc compared to the prediction. Eor B = 10^*^ G, no values of fc are excluded. Eor B = 10^'^ G, 
the cascade fraction is less than 0.18 at the 68% confidence level (CE) and less than 0.35 at the 95% 
CE. Eor B = 10^*^ G, the cascade fraction is less than 0.36 at the 68% CE and less than 0.78 at the 
95% CE. Based on predicted values and the measured upper limits for the cascade fraction, IGME 
strengths in the range (5.0-10.0) x 10^^^ G are ruled out at the 95% CE. The reader is reminded 
that this constraint only applies to the range of field strengths considered, and fields strong enough 
to produce pair halos have not been ruled out by non-detection. 

The model-dependent limits show similar sensitivity to measurements from other instruments. 
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Figure 2: Confidence intervals (68% and 95% CL) on the width of the 9^ distribution for a range of cascade 
fractions, using the source lES 0229+200 and assumed IGME strengths B=10^'^, 10^'^, and 10^'® G. 

In particular, H.E.S.S. ruled out an IGMF strength in the range (0.3-3)x 10^^^ G with slightly 
different model assumptions [8]. 

4.2 Model-independent limits 

Model-independent upper limits on the flux due to extended emission are set for all sources. 
The bulk of the beamed emission is expected to fall in the range 6^ = 0.0-0.01, thus the excess 
counts due to extended emission was calculated from the difference / — f with integration 

limits 6^ = 0.01 to 6^ = 0.24. An upper limit on the y-ray rate due to extended emission was 
calculated using the method of Rolke [14]. Translating the upper limit on the rate to a flux upper 
limit requires an assumption for the spectral index of the cascade emission. Two possible values 
were used, the assumed intrinsic index F, and a softer spectrum with spectral index (F + 1). The 
resulting 99% confidence level upper limits are shown Tab. 4. 


5. Conclusions 

A search for source extension due to cascade emission broadened by the IGMF was performed 
using long VERITAS exposures on a number of hard-spectrum blazars. No indication of extended 
emission was observed. Model-dependent limits were set on the fraction of the total emission due 
to cascade emission, assuming IGMF strengths between 10^^^ and 10^*'^ G. IGMF strengths in the 
range (5.0-10.0)x 10^'^ G were excluded at the 95% CF, assuming an IGMF coherence length of 
1 Mpc. Model-independent limits were set on the flux due to extended emission, resulting in 99% 
CF upper limits of 0.17-2.69x10^*^ cm^^ TeV^* s^* for an energy range between 160 GeV and 
1 TeV. 
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99% CL(r) 

99%CL(r+l) 




Mrk421 

0.62 

0.50 

Mrk 501 

2.69 

2.16 

lES 0229-^200 

0.64 

0.58 

H 1426-^428 

1.23 

1.18 

VERJ0521-F211 

0.20 

0.17 

lES 1218-^304 

0.58 

0.46 

PG 1553-^113 

1.83 

1.48 


Table 4: Model-independent limits on the flux from extended emission, assuming cascade emission with a 
spectral index matching the intrinsic spectral index E and (E-tl). 
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